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ABSTRACT 


The optimal control problem of a typical nuclear reactor 
power plant, which is described by a ninth-order nonlinear 
differential equation, having time-varying parameters, is 
considered. The nonlinear model complicates the optimal 
Sentroller symtnesis. Therefore, the appreach of this work 
is to approximate the response of the reactor system by 
that of a second-order linear model. The model parameters 
are chosen to minimize the derivations between the system 
and model reSpONses Ust@emaesecarchmmeutine. The optimal 
feedback parameters computed for the second-order model is 
used for suboptimal controi of the system. The model paramn- 
eters are updated to reflect the system nonlinearities as 
well as changes in the system parameters; the corresponding 
SomerouNscnemc 15 adaptive. it is shown that for the operat- 
mienGomadittons Considered, the adaptive controller need not 
be on-line. 

Also, investigation of the effects of different weighting 
HCP miMmeeleNGOStetinetion, and the effect of various control 


rod configurations on the system response are presented. 


PY ts. 





Je 


Lie 


a. 


Teper CONTENTS 


WhePROMUCTION - ="="=8- - - - - - - - - - = = 


SYSTEM REPRESENTATION- - - - - - - - - - = = 


DS i Sa 


eeOwmOh oh Niet |@omeieenee = - - = - - - - - 


ieee URE eREAGEIVity = 


i 
1 
i 
! 
! 
I 
I 
J 


ABSORBER ROD - - - - - - - - - - - - - - 
DRIVE MOTOR- - - - - - - - - - - - - - - 


OVERALL SYSTEM - - - - - - - - - - -- - 


Pree Cee tOR OF SY¥YSPEMM= = - - - - - - - 


Ne 


Bae 


oe 


IDEs 


Chet lOCt ene otomes OF SNWG@EEAR REACTOR- 


Powter CONIROL=ROD TO fis RESPONSE 
Oe Gi eee hOR = —- =9=9> - - = - = - 


Pair etek t URE COEFEICIENT- - 


Cele mrewN SOr REACTIVITY INSERTION - - - 


Crier eC GOliROUNE ROBLES = = === = = = - = = = 


A. 
se 


OPTIMUM APPROXIMATION OF NUCLEAR bmpeielees 


OMe ie CONDROL LAW= = = = = = - - = = = 


APPLICATION OF A SECOND-ORDER MODEL- 


Pie we oncCO NUS ORDERS MODEL- - - - - - = 


ee 
IB). 


eee POO tHOD= 9 - = - - = - - - - 


OPTIMUM APPROXIMATION USING PATTERN SEARCH 


So LEM 


PV AL GUESS OF eS OF A SCO 


Cetin i= = = =< - 2 os 5 
Coord ER PROGRAM OF A PATTERN SEARCH - - 
Petite ONeSOr THE RESULT - - - = - - - 


Coen LON OF A THIRD-ORDER SYSTEM- - 


10 
Re 
12 
ac 
14 
LS 
16 
19 
19 


1s 
28 
$0 
35 
oe 
41 


43 
43 
46 


47 
Syl 
a2 
D2 





VI. SUBOPTIMAL CONTROL USING A SECOND-ORDER MODEL - - 63 


WIly) REALIZATION OF ABMPTIVE CONTROLLER- - - - - - - - 73 

PeCROCEDII EP OEMINAPTIVE CONTROL - - - - - - - - 73 

Pe re TO MmenmnmemivE CONTROL - - - - - - - 78 

TUPPINIERUEN@E Or COST@EING@UION- - - -- ------ - 88 

ipo GONGLUGTONS AND RECOMMENDATIONS - - - - - - - - - 92 

PUEOUGMNeGe ween - - - - - - - - - - 92 

Pe EGG MEAT ONG ee y= - = - - = ----- 94 
APPENDIX A: NUMERICAL VALUES OF REACTOR MODEL 

ee ees eee =| - ~~ - - - - 96 

APPENDIX B: DERIVATION OF EQUATIONS ON ky AND k,- - - - 97 

APPENDIX @s OUTLINE OF PATTERN SEARCH - - - - - - - - -101 

PU BLIOCR SD Digess — 2s 2° 2 ene 

106 


INITIAL DISTRIBUTION LIST- - - - - - - - - - - - - = - 


6 





ProteOre FABLES 


COP AR SON SOLSSUMESREPERENCE POINTS 
PO VARn@wS INITIAL POWER LEVEL - - - - 


COLO ParOUeCONFIGURAGIONS - - - - - - 


PO rel) Ole: PROpeIMATED: SECOND-ORDER 


MODEL 1O VARIOUS INITIAL POWER LEVEL - 


CONTROLLER eee lS a TESTS leat 


POWER LEVEL- - - - - - 


COMPARISON OF SOME REFERENCE POINTS 
BETWEEN WITHOUT CONTROLLER AND WITH 
CONTROLLER = 5 SS 


PARAMETER ChiWGES DURING POWER LEVEL 
TRANSITIONS (50-100% FP) - - - - - - - 


-— 


Pee ERE NGio OF SMODEL AND CONTROLLER 


DURING TPGRER LEYEL TRANSITIONS (80-1003 


FP) 





TABOO SiibOLS 


Definition 


neutron or power level, as a fraction of full 
power 


time —_ 

meaCLtivity 

total delayed neutron fraction 

prompt neutron lifetime 

decay constant of qn neutron precursor 
GOneenenation Of jee neutron precursor 
fraction of delay neutron due to jt) precursor 
time constant of temperature effect 
temperature coefficient 

temperature reactivity 

reactivity due to absorber rod 

Steady state value of ékc at t=0 

time period of interest 

motor velocity 

apmplfier gain 

PliicmeOnctante Of MOtCOY 


power error 





ACKNOWLEDGEMENT 


To my advisor, Professor A. Gerba, [ wish to express 
My simecre appreciation fom the invaluable instructions and 
pidemeeniso, | give special thanks to Professor Dong H. Nguyen 


fom his advice. 





I. INTRODUCTION 


Eiecmiieeadvent Of muUcClcar reactor technology, safety 
is generally a foremost consideration in the design of a 
iol there ror and itS Gontrol system. To attain the 
PeaposceOL sate Operation With some improvement in perfor- 
Manee, there 1S a2 possible need to apply optimal control 
Pieonwe 4s COlsiderable amount of research effort has been 
iiGeectca toward the realization of optimal control systems, 
oc eteiiy tOf such big-Scale projects as the operation of 
a nuclear power reactor. It is reasonable to assume that 
ivGloateepovwen stations Of the future will be operated 
MierGetysunaer the control of disital computers, which 1s 
accepted as the most suitable device to perform such a task. 

As mentioned by Sinha and Bereznai [Ref. 1], the over- 
Simplified representation of the reactor dynamics, usually 
the one delayed neutron group model, is not satisfactory 
for control purposes, and their responses deviate consider- 
ably from that of the actual system. Therefore, the attempt 
is made to evaluate the response uSing a more complete 
representation of the reactor kinetics, and to include an 
adequate description of the controller mechanisms. 

In this paper, a realistic system of a nuclear reactor, 
which is a ninth-order nonlinear and has time-varying 


parameters, is used. The measured quantity to be controlled 





is the neutron or power level. Though any method that 
resolves the optimal control problem for this system is 
mot knewn, 1 may be feasible to get a controller that 
Will approximate the desired optimal response sufficiently 
close for practical purposes. A method that realizes 

such a suboptimal performance is described in this paper 
based on the work of Bereznai [Ref. 2]. It uses a second- 
Seader model to represent the reactor and ihe controller 
mechanism. The optimized parameters of the model are conm- 
puted to approximate the response of the actual system 
using a search routine. Then, the feedback parameters 
required to use the suboptimal feedback controller for the 
Mitweametede Orvare Casily computed. In this manner, the 
Control System is adapted to compensate for the nonlinear 


system characteristics and for the changes in the system 


parameters. 


ae 





iio ome eer RESENTATION 


The mathematical model for the plant under consideration 
momaeveltoped. To retain Sufficient accuracy, the reactor 
kimetilcs represented by the six delayed neutron group is 
ise ine Unity teedback 1S used to study the closed- 


loop response of the system. 


fee REAGTOR KINETICS 

tiewmine ttc e€QUatitons Of a chain-reacting pile have 
been derived in the literature many times [Ref. 3]. The 
Six delayed neutron group, space independent, source free 


model is described by the following equations: 


; ; 6 
anit) . eee 7(t) + 2 eerie) (Zan) 
dC. (t) Bs 
<a = oie n(t) - d. C. (t) ((2322)) 
whlemesd = fo 2, «..; 0. 


PeelpOGetnese cquations are linear in the so-called state 
Mesables, n and C. (i= 1, ..., 6), and linear in the input, 
Meee DUtethe syscem 1S not jointly linear in state and 
control. Any feedback of n to formulate ék(t) results in 


a nonlinear system. 





foe LehiPERATURE REACTIVITY 

In most reactors, heat-transfer dynamics are coupled 
Pomme OnetInetlesmby a temperature coefficient of reactiv- 
ity which is generated by thermal expansion of the core and 
by change in various neutron cross sections. In other words, 
ioe eancton reduces densit,, amd thereby the modera- 
moi bitdty and reactivity. Also, the neutron capture-to- 
fission ratio of the fuel is altered by temperature changes, 
and a smaller effect, due to changes in absorption and 
fission cross sections, increases leakage and decreases 
reactivity. Further, thermal-reactivity contribution 
results from so-called Doppier broadening associated with 
the resonance region of the neutron energy spectrum. The 


temperature coefficient a, may usually be assumed to be a 


t 
Semotant Within the accuracy with which it can be pre- 
dicted. Hence, the temperature reactivity may be expressed 


2S 


a et Ent) ; (2.3) 
t 1E 
The significance of the temperature coefficient of reactivity 
is obvious from Figure 2.1. Even if a linear heat-exchange 
model is wediizedl the temperature-reactivity feedback 
wie meeiemheadectoer dynamical model to be nonlinear. The 
effective temperature is assumed to be proportional to the 


operating power level. 





Ne S<) os Neutron n 


ranecics 


AVS yes suRtbae 


CGOeCrEriIcilent 





FAO Be zee) 


Feedback Loop cf Temperature-Reactivity 


C. ABSORBER ROD 

The reactivity term 6k is the sum of the externally 
applied reactivity due to the absorber rod, and the change 
Mute meactivity due to temperature varmations. The mathe- 
matical representation of the reactivity due to the absorber 


rod is given in Reference 2: 


Sk_ = 6k, + S'W(V) dt (2.4) 
E 0 0 
where W(V) = 0.02 V(t) igh: elas (a3) 
= 0.3 Neel 
= -0.3 WO eas 


Boat 1S: 
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POM Gemcee2} 


i 


Amseaper kod Characteristic Curve 


W(V) has the units mk/sec, and represents the reactivity 
Sramee auc tO the motion of the absorber rod. This formu- 
lation assumes that the control rod is near the center of 
mts movement, where the reactivity changes linearly with 


meeance trom the Center of the core. 


D. DRIVE MOTOR 


iemerme response Of the ™drave motor is given by: 


dV 


ratte (220) 
i 
where Cromartie 15 MOtOr velocity, 


€ iS power error and am SS eae 


comseent:t Of MOtCOT. 
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POE RALL SYSTEM 
The block diagram of the nuclear reactor model considered 


in this thesis is represented in Figure (2.3). 







Ae Billet 












MGMOn. Lx Group Power 
Decor Neutron ——— 
Rime tics 


Lemperatur 
Coefficien 


Fagure (2.3) 


Block Diagram of Nuclear Reactor and Reactivity Mechanism 


In order to study the closed-loop response of the system, 
unity feedback is used. The error signal (e) that expresses 
the deviation between the demanded and actual power levels 
momaimtpliricad, and the output iS applied to the absorber rod 
Gmive motor. The rate of reactivity insertion (Sk /sec) is 
pueportional to the error signal until the motor reaches 

its maximum speed. The actual reactivity 6k, which is a 
MedoWme uO the Criticality cr multiplication factor of the 
Teactor, is the difference between the reactivity worth of 


puemdiosonber 10d and the effect of the temperature change. 
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In state variable form the system cquations that represent 
the reactor and the absorber rod mechanism in an open loop 


configuration are written in the following: 





n ntW 
Ao Az Ag Ag Ag | BY {2 0 
, of @ @ Go we 0 C, 0 
-A, 0 0 0 0 0 0 fig 0 
ooo 8 0 C, 0 
e 
= (~ 
0 0 -A, 0 0 0 0 C, I+] 0 
0 0 O -A, 0 0 0 e 0 
0 0 0 0 -Ag 0 0 - 0 
=i ; 
oO 0 0 0 0 = 4 sk.| | 0 
T 18 
t 
0 0 0 08 0 o ct V om 
T i 
m Mh 
(2.7) 


he ok 
J 





Using matrix notation 


xo eet Be 


and the system output is given by 


where cr = i. Od 


(25. ou) 


Cao) 


The numerical values of the parameters are given in 


Appendix A. 


Hence, the parameters of reactor kinetics 


eGeattonswearc tygeical values, and the parameters of the 


absorber rod mechanism are given in Reference 4. 
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1 ipl i “DYNAMICAL Pea LOk OF SyoreM 


iiemeduatlons OL time SyStem described in Section II were 
Solved on an IBM 360, Model 67 digital computer (W. R. Church 
Computer Center at the Naval Postgraduate School) using the 
method of digital simulation language (DSL), which has been 
Suitably modified to take into account the nonlinear nature 


Grethe differential equations. 


weeeGhosmw-bOOrP RESPONSE OF NUCLEAR REACTOR 
Inewelosed-1o00p responses of the nuclear reactor to 
Step-change in demanded power level, for various initial 
power level are shown in Figure (3.]) to Figure (3.8). 
imeeeberterence pointS such as the time to reach first full 
power (100% FP), the amplitude and the time of maximum 
Svershoot are tabulated in Table (3.1). The nonlinear 


characteristic of the system is clearly shown in the figures. 


Peer Eel eOr eCONTROL-ROD TO THE RESPONSE OF NUCLEAR REACTOR 
Various control-rod configurations are tried to examine 
the effect to the system response. These are shown in 
Table (3.2). Hence, the number given EOmeieCmecOn uno) — hod 
shown in Table (3.2) are corresponding to the number of 
Emotres: (5.9) and ¥(3.10). 
iiemeomemourations of 1, 2 and 3 have the same eee 
value and different slopes for the linear region. These 


System responses are shown in Figure (3.9). In this case, 


i 
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No. 


Table (3.2) 


Control Rod Configurations 


Configuration 





ie 


W 


Equation 
DOA As LS 
05 aE 
-0.3 V < -15 
0.06vV [vi < 5 
Om os 
-0.3 V< -5 
OO |v |< 30 
68 V > 30 
56 V < -30 
0.02 |i, <7 7210 
0.4 V > 20 
-0.4 V < -20 
T.OLW |e alte 
D2 V > 10 
Ore V < -10 
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Steady State Response 
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Pitre (3.3) 


Step Response of Initial Power Level 90% 


2 4 6 8 ti 
Fionere (3.4) 
Step Response of Initial Power Level 80% 





Power 


ibe Ce 
100% SS 
80 
4 8 12 16 t 
sec 
Figure (3.5) 
Step Response of Initial Power Level 60% 
100 
90 
80 
70 
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5 0 ees eee cetera ene a Petinatarwea a er ata eater ar rp A SS 
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Figure (3.7) 
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the early system response for all three configurations are 
Me@@ledl. On the other hand, the configurations of 1, 4 and 
5 have different We values. AS understandable from 


Paoumermmes. 10), the different W give much effect to the 


max 
Seley tame response. S50, the control-rod characteristic 
becomes the important factor in determining the system 
response for the early values of time. In this paper, the 
control-rod characteristic of 1, which has an application 
Mmas particular reactor study {Ref. 2], is used in the 
mebowine Sections. In addition, a faster rod response 


Gmaracteristic 1S used in Section D where maximum reactivity 


insertion is investigated. 


Seer eG! OF TEMPERATURE COEFFICIENT 


The temperature coefficient o, is assumed to be a constant, 


16 
Diminmemd) Dracticam reactor it May vary over some range. 
To investigate the effect of some range of the tempera- 


ture coefficient, the following temperature coefficients 


are used and the responses are Shown in Figure (3.11). 


Us 


a, = -8 x oem 
Sr eae ae ie 
a, = HUge>.< 1a © 
a, = 10x 10° 


As the figures indicate, for values of the temperature 
coefficient about the nominal value of -3.5 x ‘Lees these 


temperature coefficients do not have much influence to 
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the system response. It works to increase or decrease 


the overshoot. 


DB. CONSTRAINT OF REACTIVITY INSERTION 
[opdemnederavyivy Ok which as the sum of the reactivity 
Mc rO sticeCONnErOL rod, and the temperature reactivity is 


expressed as: 


6k = 6k. + 6k (3.1) 


ieee previous Sections, the magnitude of the reactivity 
was not constrained. As shown in Section C, the magnitude 
of temperature reactivity is not very influential in the 
system response. The interested point is the constraint 
applied to the control-rod reactivity rate. The movement 
of the control-rod should be restricted to a maximum value 
from physical viewpoints. 

Now, the graph of the control-rod reactivity of the 
initial power level 80% is shown in Figure (3.12). Hence, 
gsecescribed in Section B in Chapter III, the faster 
control-rod is chosen. Three different constraints of 
Pomerol-nod ,eactivity are conSidered. These three con- 
StraintsS are interesting points, because the first works 
to reduce both the first and the second overshoots, and the 
third constraint is entirely below the actual reactivity 
except the starting area. These situations are shown in 


icine (35,12). 
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TWemereure (5.15) shows the relation of the system 
response and the control-rod reactivity. The obvious fact 
is that the control-rod reactivity gives much effect, 
especially stability, to the system response. The influence 
Of the constraints of the control-rod reactivity is shown 
Meereure (3.14). As shown in Figure (3.14), the 
constraint of OK. = Ses 10 = Homer nc mOst eLrective one, 
and the desirable system response of fast rise and minimum 
overshoot is obtained. So the suitable constraint on the 
MedctiVity wnsertion of the control rod contributes to the 
stability of the system while improving the system 


performance. 


a2 








Sia 


AYTATIIEAY 


ALTAT}OVOY pue ssuodssy weiskSs FO UOTIETOY 
Das (CL°¢) eansty 


c Z —_ 


AYTATIIVSY 


osuodsoy woysdAs 


Oe 


06 


OOT 


OOT - 08 (%) ToOAdT 29OM0g 


53 








¢ 


it x4 iG 


SJUTEIZSUOD) FO BOUSNTFUT 
el Sjeeins es 
¢ Z 


IIIS D eG 


OOT 


SHO! [3 


[28A9NT 19M0g 


- 308 


o00T 


34 





Veo reine CONTROL PROBLEM 


There is a possible need to investigate the application 
Of Optimal control theory to nuclear reactor system in 
order to improve the system performance. In this thesis 
the realization of optimal control Systems will be investi- 
gated for step demand changes in power. The procedure is, 


however, applied for any deterministic demand input. 


A. OPTIMAL CONTROL LAW 
The optimal control theory has been derived in many 
iererauures [|[Rets. 5, 6 amd 7|. The optimal control problem 


Prom ing da cCOmtero! P*eU which causes the nonlinear system 
Gye er CoGt)., W(t), t) (4.1) 


womtoliow a trajectory x*eX that minimizes the performance 


Newell € 


e 
J=h(x(te), te) + * g(x(t), w(t), t) dt (4.2) 
t 
0 


Por most Cases of practical Viewpoint, such as the control 
Meee lecarerecdetor, the Optimal control must be realized 


in a feedback manner, 


y= £(x(t), t) (4.3) 


YI 
ur 





While the mathematical solution of this problem is established, 
tminvorvyes tne SOLUTION Of @ nonlinear partial differential 
equation. 

Optimal-controller synthesis 1S complicated by the 
nonlinear model which is utilized to approximate reactor 
dynamics. For linear systems, however, it is possible 
PemoOotaimeanalytical Solutions to such problems. Further, 
PeioeonOmMe times Section that the synthesis of the control 
which minimizes an integrated quadratic in state and control 
Only requires a linear feedback of state variables. 

Assume for convenience that the dynamical behavior of 


the system is approximated by the linear model 
ee Ce) (t) + B(e) v(t) (4.4) 
where x(t) = xX 


~ 0 


The problem is to find the admissible u(t), or, preferably, 
Meoreiwiteheecentrols equation (4.4) with respect to some 


reference x*(t) so as to minimize 


£j(x(t), u(t), t)dt (4.5) 


where 


£ (x(t), u(t), t) = 5 [(x*-x)" + Q(t) (X*-x) 


tut) R(t) yl (4.6) 
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and R(t) 1s a real Symmetric positive-definite matrix and 
Q(t) is a positive-semidefinite matrix. 

Obviously, this problem may be solved by the maximum 
Poarncaple mon Hamilton's GCquations), or by dynamic program- 
ming (or the Hamilton-Jacobi-Bellman equations). To use 
the Hamilton-Jacobi-Bellman equation, the Hamiltonian is 


formed: 


H(x(t), u(t), $25 t) 


~~ 


= 5 [x! (t) CAG) Cees a" (t) R(t) u(t) ] 


as! 
+ 22°. fact) x(t) + BCt) u(£)) (4.7) 


m necesSary condition for w(t) to minimize H 1s that 


=— Oe thus 


0H oJ 
ou (x(t), u(t), ipa 5) t) 


~~ 


T J 
= R(t) p(t) + B(t) $= 0 (4.8) 
snee ie ma tt 1x 
2 
sy = R(t) (4.9) 
UW 


eeoO meivendetinite and H 1s a quadratic form in wu. the 


~ 


control that satisfies equation (4.8) does minimize H. 
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Solving equation (4.8) for u*(t) gives 


=i i JJ 
u*(t) = -R'(t) B(t) 55 (4.10) 
which when substituted in equation (4.7) yields 
i 
OJ a eso ole) =I ee eid) 
eee I) Se tox EEF x 
) 
OJ 
the Hamilton-Jacobi-Bellman equation is 
ere x , t) oot x a2 Cet) 
___ i eres a Qx- Lg eee. a B roi pi ne 
dt i 2 “ers ~~ ~ 0X 
aj See 
ee 9) 
. ee A X (4,12) 
From equation (4.5) the boundary condition is 
J) (x(t ¢) J te) = 0 (AS) 


Simce the Minimum cost for the discrete linear regulator 
PeowvletheiS 42 Quadratic function of the state, it seems 


reasonable to guess as a solution of the form [Ref. 5] 


ii 


J (R(t), t) = Fx (t) K(t) x(t) (4.14) 


where K(t) is a real symmetric positive-definite matrix 
@idteis to be determined. Substituting this assumed solution 


meeequation (4.12) yields the result 
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BR Deemer Keiox = 0 (4.15) 


where generally all these matrices and vectors are functions 
femerme.  Equatron (4.15) must be valid for all x(t). 


ism tne Matrix property, this equation becomes: 
K(t) # K(t) A(t) + A*(t) K(t) 
-] T = 
-K(t) B(t) R(t) Bo(t) K(t) + Q(t) = 0 (4.16) 


and the boundary condition is 


K(t¢) = (4.17) 


Once K(t) has been determined, the optimal control law is 


given by 


w*(t) = -R*(t) Bi(t) K(t) x(t) (4.18) 


~e 


MetnessysStem is completely controllable, it is sometimes 
desirable to let te * ™, SO that wthe optimal control 
accurately maintains the desired terminal state once it is 
meaencdea ss further, Kalman {Reft. 8] has shown that if 


oO 


D= FS (x(t) Q x(t) + w(t) R u(t)jat (4.19) 
Pee y 


where Q and R are positive-definite, constant, symmetrical 


~ 


matrices and the system [equation (4.4)] is time invariant 


So 








(A and B are constant matrices), then 


lim dK(t) _ 


eae uy 





In this case, equation (4.16) becomes 


ah 


BK + Q = 0 


neh Kk BR 


rometas problem, the optimal control is 


and the synthesis requires a constant linear feedback 


the state variables as shown in Figure (4.1). 


< (t) | x(e) 


7 at: 
6 f dt 





Controller 


Pieire (4. 1) 


Optimal Feedback Controller for Linear Regulator 
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(4.20) 


(CArce7 15) 


(4.22) 


one Bisll 





poe te ON OF A SE@OND-ORDER M@DEL 


icc cou Ondcpmametrtenmential cquation 1s expressed as 





2 
dex dx _ 


The elements of K are determined by solving equation (4.21). 
The corresponding expressions for the optimal feedback 
parameters are obtained by using equation (4.22), and 


developed in Appendix B with the result. 


a q = 
ee 2 fear 
ky = By i: Bo [an + — bo] (4.25) 
a 6 | 1 
a ae 2 D 
ky B + B [ay + = bp + 2boko] (4.26) 
where 
fy Kya 
ea , Kaa Ko= 
0 wy King Ka 
b.k b.k 
2 te _— ME 
and Ko = - ; kK, se 





4] 





and selection of the positive parameters corresponds to the 


positive definite requirement on k. 
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VY. OPTIMUMPAPPROXIMATION OF NUCLEAR 


Peano eo) salem pi A SoECOND~ORDER MODEL 


The nonlinear high-order system described in Chapter II, 


ween 1S Written below — 


x= Ax+ Be (5.1) 
where x(t) =X» 

e —-- XxX 
and eS) = cr Xx C52) 


may sometimes be approximated with sufficient accuracy 

by a linear model for small variations in state and control. 
ICmarplyetitesoptimal control@law to the nuclear reactor 
System, this optimum approximation problem is the most 
necessary thing. This idea has been developed by Sinha and 
Pomeanameykher. 9). But in this paper, the different 


Seron wcorheerla are used. 


eo VesteiianeON O&eMETHOD 
Considering a discrete set of values of men taken over 


a suitable interval of time, 





Y = ty}, a>? pene s y.} (S.3) 


This set represents discrete sample points of the response 
Ou the system, which is obtained by solving equation (5.1), 
using the digital computer. These results are shown in 
Captor ity The continuous output response Yt) 1s sampled 
at sufficiently close intervals of time lest the significant 
information is lost. The objective is to find another 
output set y* in optimal fashion, associated with a second- 


order model described by the equations. 


eae x BS a (on. 4)) 


~ ~~ ~~ 


More Ca (5.5) 


such that, for the same input, the following objective is 
Satisfied. 


A scalar error performance function J is minimized 
J = g[W Cy, - y#)] (5.6) 
edie PT: 1 ; 


which is some suitable function of the errors (y; - y=) with 
a vector weighting factor W. attached at each sampling 
instant. 

The choice of the error criterion expressed by 


equation (5.6), has a direct effect on the parameters of 
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the approximating model. Since the purpose of the error 
Suite nloneis sto Medsure the extent to which the response of 
the model deviates from that of the actual system, the main 
Dpaeebeemr iS Gow tO express this déviation numerically. It 
is normal practice to consider the norm of the output error 
paaeocaseeit to Some power p. This objective becomes to 


minimize a2 summation of deviation. 


I2 (S273) 


for the single-output case, where {fy.} is a scalar sequence, 
1 

ama the weighting factor We 1s unity, equation (5.7) means 

a measure of the area between the curves when P=1, and the 


mean square error when P=2. 


Response 


alll a ee me 
rd big 
~~ a VY 
ies. | > ~---~ 
t 
Sa yt (th) 
/ 
/ 
= 
aN 
ih hs A Ae) 


haotunre (S21) 
Reference Response of the System 


and the Optimal Second-Order Model 
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Consider the system response a) and the model response 
te ae NovwmMernerioutre (55h). Et 1S required to find a 
model of the system so that the error criterion iS mini- 
mized. The following two error criteria as the performance 
measure are used to approximate the response Y(t) (WO 18 als, 
response Vener of some model. 

fh) Minimization of 


Z 
= ~ yk 
Jy : OF 2) (5.8) 


(Aaidnimization of 


See - yt)? + (ayia. a74)2y (5.9) 
7, OS ade dt dt 
by considering the derivative of Va and YG: the more signifi- 


cant measure 1S attained where the slope is changing rapidly. 


B. OPTIMUM APPROXIMATION USING PATTERN SEARCH 

The problem of approximating high order system by a 
second-order model in an optimum manner can only be solved 
Pope elite Cmnor Griteria by the presently available 
techniques. 

Since an analytical solution to this problem does not 
appear feasible, various search techniques may be considered. 
While gradient methods are generally quite efficient in 
Pomc wanilininlies the Necessity Of finding the partial 


detivatives Of varbitrary error criteria with respect to all 
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the model parameters, becomes a great disadvantage. 
Pimeiewounem Nand, direct search techniques involve 
evaluating the effect of sequential parameter changes in 

an organized manner. The pattern-search technique of Hooke 
and Jeeves [Refs. 10 and 11] was selected as a suitable 
Metiod ton this task, Appendix C shows an outline of the 
pattern-search technique, which is contained within the 


IBM 360/67 library under the subroutine name of DIRECT. 


fee er AL GUESS OF COEFFICIENTS OF A SECOND-ORDER MODEL 
An important problem associated with search technique 
is the selection of the starting parameters, since these 
have considerable influence on the convergence of the 
process, and on the probability of locating a local 
Minimum of a performance measure. Therefore, the reasonable 
initial guess to compute the parameters using subroutine 
MmiReGiwaseneeaecd. ithe characteristics of the second-order 
linear system have been developed in many manners [Ref. 12]. 
Since the response of a nuclear reactor system 1S approxi- 
Maced by a Second-order model, the characteristic of that 
mS applicable to determine the coefficients. The differen- 
ibtateeaquation of a Second-order system in the form of the 
mative! £requency ,; Wa and the damping ratio, 7, is expressed 


as 


KU(t) GG) 





ee) 
oe 
Dh 
SX 
= -4 
eal 
os 
— 
—, 
va 
tt 


KA 
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The response x(t) and its first derivative are obtained by 


standard methods, and are shown below. 


ee re Gin (Wl - 7c? t-w) cole 
a: ——— n 
n Late 
-cWt 
a) = See 71 - ¢* ¢) eS 12) 


mene A o21S see plet of x(t) and mi Since the aU 
1s a unit step, the value of A equals to one. To make the 


Steady state value of x(t) unity, 


x(t) = Sy = ys Cts) 


Therefore, if K = W_* 


, the steady state value of x(t) 
becomes unity. 

The value of Xnax may be found by a setting equation 
(le) to zenro, finding ty SOMSUNS telecine this value of 


ty imto eqwet vem (5.11) 


-cW t ean. 
x(t.) = 0 = ats eB sin (W, ES (5.14) 
18) ee r 
Then 
we = TI CS)..15) 
my Pp ; 
Cee 


4 













Preune (5.2) 


x{t) and x(t) vs. time 


and 


= = G 
ee ee ot CS 26) 


for a special case of a second-order system expressed in 


Sauderon (4.235), 


Xo = ag X, - ay X, + by Le) 


a 7 2 = 
As stated previously, choosing ay = bo Ne » ay 2cW 


and U(t) = unit step input, the steady state value of xy 1S 
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oO 


a(S ee aS esas) 


0 
The characteristic response of the reactor system of the 
initial power level 90% solved by digital computer is shown 


tieereure (5.3) 


Power 
Level 


(6) 





Hac usemmeo.. 5) 


Response of Nuclear Reactor System (90-100%) 


Subs tieu ting t, = 2.48 and ewes > iO ean co Cala ions 


(5.15) and (5.16), the computed natural frequency and 


dampine ratio are: 


—— 

<, 
JX 
it li 
Nh (eae 
-—_ ee) 
NO eae 
OO oe 
wn 


Then, ag = W = 4.528 
n 


0b) 

I 
tO 
ae 
a 

It 
WGN 
ix 
ho 
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Using these values as the initial guess coefficients, the 


computed optimal parameters by search routine are: 


ier 54 


Q9 
i! 


ay feral 


there 1S considerable difference between the initial guess 
values and the computed optimal coefficients. In a practi- 
Cal point and considering the nonlinear nature of the systen, 
these initial guess values are sufficiently near the 
PieimeoctrereitUsmmoimce there 19 a definite necessity 

to compute initial guess values, this method is shown to 


be acceptable. 


Pee COMPUTERS PROGRAM OF A PATTERN SEARCH 

A computer program has been written that uses a pattern 
search subroutine to find optimum second-order models for 
high-order systems. The program has the following summarized 
features.. 


fjminesougemG responses of the nuclear reactor system to 
unit step function are available at discrete uniform 
intervals of time. 


(2) The parameters of the model to be used for starting 
values are given by the user. 


(3) The program assumes a uniform weighting sequence. 
(4) The second-order model is solved by the Runge-Kutta- 
CtemeoMeri-Ooncdcer metmod, Ail] calculations are in 


double precision. 


(5) The objective functions to be minimized are 
eeuaerons (5.6) and (5.9). 
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E. MANIPULATION OF THE RESULT 

DconomipiemsSearcen rOUtane, the coefficients of the 
approximated second-order model are manipulated with the 
results shown in Table (5.1), and the plot of that is 
Eumoure (5.4). 

The difference in the model coefficients using the 
performance measure Jy and that of the performance measure 
J, fool tee lines 15 elearly shown in Figure (5.5). 

The system responses and the approximated second-order 


model of the performance measure J to various initial power 


levels are shown in Figure (5.6) through Figure (5.11). 


Eee CONSIDERATION OF A THIRD-ORDER SYSTEM 

Comparing the system response with the second-order 
model as shown in the previous part, there is some appre- 
ciable difference. But, as the considered system has a 
ninth-order nonlinearity, to approximate that by a second- 
order model, the difficulty exists in the natural sense. 

Pomtlacmolrrcnrence Is CONSidered acceptable. As a trial, 
even though the optimal control law of the third-order model 
1s not proposed, the approximated third-order model is com- 


puted using the same search routine. The equation of that 





Tee 
3 Z 
Geax case Gite = 
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Taple (5.1) 


Coefficients of Approximated Second-Order Model 


to Various Initial Power Level 


Initial Power Level Time Interval 
Nic = 0.95 (28 sample points) 25 
Jy = 8.191 x 10° 
ao = 4,399 
ay = 2.149 
JB Me IO se ie, 
ay = 4,008 
ay = 1.809 
Nic = 0.9 (28 sample points) OS 
J, = 6.60 x ron 
ay = eos 
ay = 1.614 
= ~3 
J, = Seed x. 10 
ay = 1.554 
ay = 1.569 
Nic 7 O38 (28 sample points) O25 
5, = 3.792 x io 
ay = 0.507 
a, = 0.950 
Jo = eee 10° 
a = 0.478 
ay = 0.885 
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Lc» 


ne 


ne 


(28 sample points) 


1.791 x 10 7 

a = ec 

ay = 0251) 
-2 

3.529 x 10 

ay = 0.137 

a, = 0.496 


(28 sample points) 


noes 0 
a = 0.0804 
ay = 0.358 
-2 
6405. x 10 
ag = 0.0795 
ay = 0.3551 


(28 sample points) 


Oe enions - 
ae 2 elo x Lg 
go) EGE) ox LG! 
1.409 x 10° 
Diy SMa MIO XS NO 
eo Sx 10 
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(28 sample points) 


6.097 x 10% 
ay = 1.564 x 10° 
a, = 6.097 x 10° 
3.560 x 10°} 
a SEL on 
a, = 6.966 x 10° 
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Promien( 5.4 ) 
Ccomaretents Plot of Second Order Model 
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99 
~~ System 
98 
97 
96 fp 
95% | | 
l d 2 : by sec 
Figure (5.6) 
Comparison of Step Response of System 
and Model to 95% Initial Power Level 
Model 
_ Lea 
"Sys tem 
90 
80% 
2 4 6 8 
Sec 


Eaiteume. (5.7 } 


Comparison of Step Response of System 
and Model to Initial Power Level 803 
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Model 


ae 


100 an = ao es 
8 0 = ae 
System 
60% | | 
4 8 eZ 16 t 
sec 
Figure (5.8) 
Comparison of Step Response of System 
and Model to Initial Power Level 60% 
— > 
100 (a a5 
90 
8 0 
feaz———-—- System 
70 
60 
50% a. 
4 8 az 16 4 
sec 


Fissure (5.9) 


Comparison of Step Response of System 
and Model to Initial Power Level 50% 
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Oe 








Figure (5.10) 


Response to Initial Power 30% 
100 —— : » SS 


AS 
80 N, 
Model 

60 

oystem 
40 
A | 

4 & FZ 1G 


Figure (5.11) 
Comparison of Step Response of System 
and Moacl to Initial Power Level 20% 


60 


20 





where U(t) = unit step input. The obtained data of the 


optimal third-order system to initial power level 50% is: 


Jy = 6.213 x ‘V0 
ap = 25 10 
ay = 8.10 
a. = ise 
a, = 10.64 

> 
v2 oye. 
Z 

ay = 2 Ao 
a) = 8.53 
ao = 10.63 
az = 10.88 


iMemelarittedsplot 1S showm in Fissure (5.12). Also, the 
second-order model is added in this graph. It is observed 
wicdiemene —Ctilira-opdeq = model iS more realistic than the 
second-order model to approximate the actual system. Since 
temas desirable tc construct the optimal control law of 

the third-order model, that work will be expected in the 


MELT; < 
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VI. SUBOPTIMAL CONTROL USING 


A SECOND ORDER MODEL 


The suboptimal control of a reactor system using optimal 
second-order model is made possible by the fact that the 
State variables available for feedback are the actuai output 
power level and its rate of change. The optimum linear 
feedback of these variables for the case of the second- 
Order model with an integral quadratic cost function exists, 
and the appropriate feedback coefficients can be determined 
explicitly in terms of the parameters Ag» ay and the 
constants used in the cost function. These realize a sub- 
optimal control for the reactor system. The block diagram 
representation of the suboptimal controller, the reactor 
System and the optimal controller for the second-order model 
sesshowm im bPieures (6.1) and (6.2). 


same taple (5.1), the calculated k, and ky are shown 


0 
Pipmiamhow( Oe jie Hence, the weirghting factor q, = 4,77 in 
equations (4.25) and (4.26) is used to compute kg and ky. 

In this system and model with the suboptimal controller 
and the optimal controller, one modification is needed 


to make the steady state value unity. 


Prem Paeure (6.2) 


xa (Ss) a 
sy 7 > (G25) 
ne Se hay aS 7 Ay 
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habe (oe) 
Cont rower soonmaneters 


tomycserous lwreaal Power Lével 


Initial Power Level k k 


0 JI 
95% 0.4142 02736 
90% Oy aakaZ 0.602 
80% 0.4142 OF 6.51 
60% Dea 2 0.869 
50% eS, Tegel eal 
and 
U(s) = k U, Cs) - (ky ie kjs) x, (s) (on) 


PMbStmeucino (oO,2) into (6.1) 


1 
x, (Cs) ak 


go Ono) 
uy (s) 5? + (a, + kj) S t+ ay Ol kp) 


Using the final value theorem 


_ lim : 
Se) 7S semey. x, Cs) 
_ 2. yk 
2 9 1s 0 (ford) 


DMetererc., the t¢eed forward controller of the amplitude 


k=) ] + Ko to make the steady-state value unity is needed. 
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To make the comparison possible, between the responses 
Spee sumoOmtimaleconcroller and the corresponding optimal 
model controller, and that of the system and the model 
Teeotimany scomurollen, Figpure (6.3) shows a plot of the 
System response to initial power level 90% without the 
Sontrollcreand with the comtroller. The describable fact 
is that the maxiinum overshoot of the system response with 
MicmoInoOvEmmalmeontroliler 1s reduced. Also, Figure (6.4) 
Shows 2 plot of the model response to the same initial power 
Tove itmwe etme Controller and with the controller. The 
desirable reduction of the time to reach the 100% power 
Revelmiswepservea in the model response with the optimal 
Senor lon. sommmtne Considered suboptimal controller 
works to make the reduction of the maximum overshoot of the 
system response with no effect on the early time response. 
Cette OGhcrenane. the Optimal controller of the model 
does affect the early time response by making it faster. 
Table (6.2) is the comparison of the overshoot and the time 
to reach the first full power (100%) of the system and the 
model, without controlier and with controller. 

The responses between the system and the model to 
various initial power level with the controllers as shown 
Mitiremom 7 lethrouen Figure (6.9). Comparing these 
(Parienettem( Sosa) through Ficure (5.9), the controllers work 
to effect the model rise time in such a manner to make the 
overshoot of the system and the model occur at the same. time, 


and have approximately the same valuc. 
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Without 
102 Lo SGinteom lea, 
ae 
100 LO Sea 


4 aon 


98 
With 
Controller 

96 

94 

92 

90% 

dt Z S 4 5 
Pacire (655) 
Comparison of System Response of Initial 
Fovereleverl 904 without Controller 
and with Controller 
Without 
Controller 
100 Fee ee 
98 
eee WN Teal 
Comecoller 

96 

94 

92 

90% | 

1 zZ 3 4 S 


ERigune™ (6.4) 
Comparison of Model Response of Initial 
Power Level 906 without Controller 
anaewe tin Controller 


69 





Model 


100 ——— 
99 
OR 

ee —~ System 
~ 98 
ro 
= 
o 97 
c 
© 
= 96 
O 
Ay 

95 ) 

it 2 g 4 t 


Sec 
Figure (6.5) 


Step Responses of Svstem and Model 
with Controllers to Initial Power Level 95% 


Model 


100 


aS 


Power Level (%) 





90 


eo FE EE a ne Se 
1 2 5 4 


Cec 


Figure (6.6) 
Step Response of System and Model 
WiecentnOolrlers to Initial Power Level 90% 
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Jee ke ZAal ton Or ADAPTIVE CONTROLLER 


The effect of the suboptimal reactor controller and the 
optimal model controller is described in the previous 
nineteen c une COnurdm@ler parameters are fixed. It 
is possible to design programmed time variations of 
SOC cOuben parameters tO achieve instantaneous near optimum 
Pero wEOnetnessyStem at alll times. In addition, the 
Soremepagamerzer Such aS the temperature coefficient 
change with time and the existence of such conditions is 
the reason for application of an adaptive control system. 
Eveleigh [Ref. 13] has defined the adaptive control as 
follows: 
"If an index of performance (IP) 1s available which 
indicates the system's instantaneous or short-term 
average performance quality, and if a control loop 
is set up to optimize the IP automatically by adjusting 
controller parameters, the parameter-adjustment configura- 
Elon iswecalled an adaptive control loop. It is important 
Pov wWidmedieadaprine controller is a parameter- 
adjustment loop above and beyond the normal feedback 
Mesa NrOrcOnenolmpostevon, velocity, and the like. 
Neaperve scOmerol 1S thus an effort to extend basic 
Peni COMenel "cOoncepcsS tO time-varying Systems." 

A. PROCEDURE OF ADAPTIVE CONTROL 

Hetlecwetmeweconeept Ot the adaptive control, which is 
shown in Figure (7.1) as used in Reference 2, is employed. 
This procedure is described below. 

The model parameters are precomputed, over the expected 


range of demanded power level changes, to be uSed as starting 


values for the optimization. As the model coefficients 


io 
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rarounre (7.1) 
Block Diagram or Adaptive Control System 
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are updated, the parameters of the feedback controller are 
computed from equations (4.25) and (4.26), depending on the 
(iirc Ountine COSE tUNneLtlon, = Here, also, the weighting 

fae cor q, = 4, = 7 is uSed to make the comparison possible 
MPiemmineirols Cuanter,  Ihat 1S, the cost function has the 


moO lowime tOrm: 


CO 


cae a Ce 
25 ~ 


>» Uo) 6B (7.1) 


Initially, the nuclear reactor is assumed to be operating 
at a steady state power level. At time t = 0, a step 
ei2ngeurO a Newm operating Jeve] 2s required. Using the 
model parameters appropriate for the demanded change and 
PaewaesitreadscOst 1unction, the controller parameters are 
evaluated, and the system begins the transition to the new 
power level. Over several sampling intervals, the value 
of the system response at each Sample point is stored. 
This means neither updating of the model nor recomputation 
of the controller parameters takes place. 

The fitness between system and model depends on the 
a@etracy Of the Sterting model. The model 1s reoptimized on 
the basis of the observed system response. The next step 
1S to recompute the controller parameters. After completing 
Pues seepmoccaures, tie most Optimized controller parameters 
PeecoOlcmimecr Val 1S f£Olmd, This optimized controller 
introduces the desirable system response. Repeating these 
steps, the system response of entire intervals is obtained 


im the most optimal manner. 


—_— 





Asedesceriped in Section B of Chapter III, the high 
gain of the control rod rate gives the response time 
Hieiminkiwence. lo make faster response, the following 


eomenolerod 15 USea instead of equation (2.5). 





emer Cre te/ 22.) 


Meoheoain Gontrol Rod Characteristic 


The system response without controller of initial 
power level 50% using the high gain control rod is shown 
in Figure (7.3). Also, the approximated second order 
model is shown in the same plot. The time to reach 100% 
Malue Or the system 15 eCarliemewith the higher gain rod. 

Apparently, judging from this graph, the second-order 
model does not fit to the system response. Though the 
third-order model describes the system more accurately 
than the second-order model as shown in Figure (5.12), 


the usage of the third-order model is not established 
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miceHeetedlhy ana therefore iS not applied to this 
problem. 

eure (7.4) shows a plot of the system response with 
the suboptimal controller to initial power level 50% and 
Pid etouL Contmrolier to the Same initial power level. 
The suboptimal controller based on the optimal control 
law of the second-order model has the great wae to 
reduce the overshoot and the system response reaches 
Steady-state value quickly. Apparently, the uSage of the 
linear feedback controller is recommended to control the 
nuclear reactor to reduce the large overshoot and to 
reach the steady state faster. Figure (7.5) is a graph 
of the responses of the svstem and the model with the 
Contnolliens  sGemearine this with Figure (7.3), the 
second-order model with optimal controller more closely 
represents the reactor system than when no controller is 


used. 


Poe ono Oe OnmeAnAPTI VE CONTROL 

As described in the previous section, the procedures 
Gt the adaptive controller are applied by using an obser- 
vation interval of 15 seconds duration and the adaption 
intervals of 0.25 second. The obtained parameters 
Midieinemeamputcd repeatedly until sufficient accuracy is 
attained, are shown in Table (7.1). 

The corresponding responses of system and model with 


the adaptive controller are shown in Figure (7.6). Also, 
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risumemvea7lectows the System responses with the non- 
aodamenve controller and with the adaptive controller. A 
comparison of Figure (7.5) and Figure (7.6) shows that 
there 1S some improvement between the model response 
each the optimal controller and that with the adaptive 
Comeuotlero win Fieure (7.7), there 1S no appreciable 
difference between the system response with the non- 
adaptive controller and that with the adaptive controller. 
This means the concept of adaptive controller is not 
effective as long as the system parameters are constant 
and unity weighting factors are used in the cost function. 

As a reference, the response of the second-order model 
used to develop the adaptive controller is shown in 
Figure (7.8) with the system response using the adaptive 
comuroller. 

hOomoionme= tien creect Of adaptive controller 
clearly, another example is studied. The theory of 
adaptive control is applied to the system of initial power 
level 80% shown in Figure (7.9). The tabulated result 
mseohown in lable: (7.2). 

Figure (7.9) is a plot of the system responses with 
the non-adaptive controller, and with the adaptive con- 
troller. Though there is very little difference between 


them, they are almost considered identical. 
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Vii euewcr OF COST FUNCTION 


As having been described, the cost function is: 


ce 


jee pe! Oo x © Se (8.1) 
i. ae 


bo] 


In the previous chapter, the controller parameters were 
computed for the case that Q is an identity matrix and r 
equals to unity. 

Imvthis ehapter, the various weighting factors for 
the cost functions are considered, depending on the optimal 
control law. The selection of apprcpriate values for Q 
Mier rs elmp@mtant, and COcs to the core of the optimal 
regulator problem. The process of selection consists of 
Boctomines Certain values to Q and r, incorporating the 
resulting controller into the system. 

EiictiOnsoml40e sand swesZ20) O01 the controller parameters 


ares renwrittcem Lor Convenience 


i 
a q = 
_ fe a 2 1.2.2 
Ky = 5 by [aq = Pod 
A q i 
Peal &: w 2 Dre? 2 
ky = BD * BD [ay + == bo * 2bo kG] 


where qy and q, are the components of a matrix q, r is the 
relative weighting factor. 
These Weighting factors of the cost function decide 


thew controller parameters. 
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The following values are tried to check the behavior 


of the actual system. 


a 
Ca. ~ 90.1 


iiewecomemecdscONtroller parameters are: 


Ko = 0.0488 
ky = 0.1426 
4, _ 42 
oa aa * 
Ko = 0.4142 
ky = 0.797 
q q 
eee 
5) oo = 10 
Ky one raey 
ky = 3.008 
q q 
1 a Bee 
Ky ooo 
ky = 4.558 
These responses are shown in Figure (8.1). The stability 


of the system with the suboptimal controller depends on 


weighting factors of the cost function apparently. Since 
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the optimal control law is based on the quadratic cost 
function, it is essential to select the most reasonable 
cost function, as described in Reference 14. The 


resulting fact 1S consistent with Reference 14. 
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IX. CONCLUSIONS AND RECOMMENDATIONS 


The following conclusions and recommendations are 
offered as a result of the investigation contained in 


this report. 


A. CONCLUSIONS 

1. The simulation of the dynamic behavior of the 
nuclear reactor system on the digital computer to make 
feasible this investigation was obtained with some 
get fae | ty. 

To obtain the stabilized system response required 
considerable effort, because the system is nonlinear and 
described by a ninth-order differential equation. The 
unknown factors such as tne specified values of the 
nuclear reactor parameters complicated the difficulty of 
implementation of the system equations. 

een tatomeciocr. sene Saan Of the control rod, the 
reactivity-rate constraint and the effect of the tempera- 
ture coefficient have been studied. The high gain of the 
control rod causes the faster response, and the reactivity- 
rate constraint suitably chosen makes possible the more 
dese ilen LespolicesGegardine the stability, particularly 
to reduce the deviation from the desired power level. The 
magnitude of the temperature coefficient has a lesser 


effect on the system response. 





3. It has been shown that the results from the 
Opemial second-order moder of the reactor may be used 
Semecervely stor Suboptimal control of the reactor system. 
iiemsuboptimal controller of the reactor reduces the 
overshoots of the system response. Considering these 
overshoots are not desirable practically, its usage 
Pugvades an improvement in System performance. To find 
the optimal second-order model, the computer program of 
Search technique is applied to obtain the second-order 
model parameters. 

4. An adaptive control system, which may be realized 
in practice to control the power level changes of a 
iMielear redctOn, nas been proposed. 

The parameters of the second-order linear model 
are continuously updated so that the model accurately 
represents the behavior of the system. Also, the parameters 
Of the controller are updated to achieve the concept of 
the adaptive control. The result of this study shows no 
improvement of the system response using the adaptive 
controller for the operating condition considered. The 
adaptive control, therefore, would not have to be on- 
line, and could be used to update the model parameters 
to account fon Plow cihamoaio parameters Of the reactor 
system. 

Tlomvicm@ZNGinomeaetOr Of the cCOSt function gives 


the stability of the system response much influence. 
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Based on physical considerations, the weighting factor 
should be chosen as correctly as possible to provide the 
desired response and this would be a trial and error 


PHOecedure . 


B. RECOMMENDATIONS 

PeaenenevOuKerS In tis area should consider the 
following facets of nuclear reactor systems. 

Mee ticmmuctear reaeror System Studied in this report 
Gicmener cONolder the Noise problem. For the purpose of 
the present paper, it was assumed that an instantaneous 
noise-free measure of the reactor power level is available. 
However, as the temperature channel signal indicates the 
power level of the reactor at an earlier instant, becauSe 
of the finite transport time between the reactor core 
and the temperature transducer, this signal has a con- 
siderable noise component due to the turbulent coolant 
flow. 

Therefore, the noise consideration is strongly 
recommended for future studies. 

2. Since the proposed scheme relies on the use of 
Setuem rOUtIncS. Improvement in the efficiency of the 
leeie@epantieubanly as applied to digital process computers, 
would greatly enhance the usefulness of the new technique. 

3. In this report, the effectiveness of the adaptive 
control has not been shown. Consideration should be 


given to the adaptive controller for different weighting 
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factors in the cost function to find out whether an on-line 
adaptive scheme would be required. 

4. ose much of optimal control, the natural formulation 
is by state variables. A well-known and popular result is 
Pit mioreticw linear negulator problem with quadratic 
performance index as applied in this study. It has been 
Shown [Ref. 15], however, that the solution is far from 
realistic or optimum in an engineering sense. Nevertheless, 
the results and scheme obtained for the linear quadratic- 
index regulator problem have been forced indiscriminately 
on the linear time-invariant problem. 

Since the weakness of the state-variable formula- 
tion in coping with linear feedback control system design 
1s especially apparent in the sensitivity problem, which 
1s one of the primary reasons for the use of feedback, the 
sensitivity theory by Horowitz and Shaked [Ref. 16] should 


be investigated. 
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APPENDIX A 


NUMERICAL VALUES OF THE REACTOR MODEL PARAMETERS 


B = 6.41 x 10> 

B, = 2.70 x ome 

6. = 7.40 x 10°" 

Be = 2.53 x 10” 

By SL Bae ag ° 

6B. = 1.40 x OT 
Pree 10) © 

Ay = So ees sec + 
eek. 13 SGC 4 

Az = 0.301 sec 

Mp = (0. Mil sec} 

A = 3.05 x 10 2 see + 
he = )1,24>x 10 2 see t 
2 sim © see 

6 )= 950 

i: O27 Gees ec 

ee ii? 5 SCC 

ees 8,5 50 
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APPENDIX B 


DERIVATION OF EQUATIONS ON ky AND k, 
The equation is: 
Che Teer re Boeke oO = 0 
where 
0 i 0 
AR ina 
Ba) aor bo 
q 0 k 
veh ee = 11 i 
0 44 Koy = Ka 
Wal Tee Ky = Koy 
Conpucine term by term 
“ayky, Kia 7 41 %12 
Levi = 
“agk,, Koa ~ 24%22 
r “agks, “agk,, 
wage 
Kia 7 928210 a2 > 41K 22 
0 0 
caree be = 
0 
0 ae he 


.O 
~J 





bZ Kes Kk Keak 


1) 1222 

Lek baer =o 

~~ om™NlUum® ae ~ i Z 
Ko0Ko4 Ko2 


After substituting these matrices into the equation, 


rearranging, the following four equations are obtained: 


ae 

ee wo x 0 
“Ag King - 8g Kya - 7 Ky * 

b+, 
~4g Koo * Kyy - 8y Kya - 7 Ki2 Ka2 * 
b*, 
Kia 7 91 Ko1 - 20 Ko2 > 7 Kaz *22 * 9 
b*, 2 

Kig > 41 Ko2 * Koy > 8, Kan - 7 Kon + 92 * 











b2 
G3 2 
eae + = = 
medion 2g “4p 4 CP 
SO, 
-a 
0 7a ai q 
k mm ea 0 aes a2 a el b? 
eZ ) 9 y 0 r 0 
b b b 
0 0 O 


OR 





Also, 


2 
Dae 
r 
ee a 5 
b 
a) 
i 
[Mone ticwOnelmal control law, 
Ut = - RO BOK X= -k’ x 
then, 
k 
11 
1 pcs aurea | 
ko = [k, kJ oe bk > [0 bol 
ai 
b 
O ry 
1 [Ky 7! 
90, 
ee oneal 
0 ag 
ane oF 
Hl 18 





3) 








Therefore 
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REND? XC 


CUMENEs Or tii ePArTTERN SEARCH 


As devised by Hooke and Jeeves [Ref. 10], pattern 
search is a direct search method which results in relatively 
Pmteteme search alone Straight ridges or ravines. It is 
attempted with pattern search to establish the pattern 
of successful search points in the immediate past from 
Wimeiwuplausmp le LUubure Search points are predicted. The 
method for the two-dimensional case is illustrated in 
Figure (C.1). The search progresses from an initial base 
point x®, and a minimum of f(x) is to be found without the 
use of derivatives of te) A small displacement d, 
from x° in the x4 Gieection 1S ertrected, and £(x) + d,> x5) 
is evaluated and compared with E(x); x5). If the latter 
is smaller than the former, a small displacement in the 
-Xy direction is effected, and f(x) - dy; x) is compared 
with f(x), x5). Pemticelaeues 1s Greater than the former, 


a small displacement d. is made from (x4 - dy, x5) in the 


Orne c TlOn eam (x - dX ) 1S compared with f(x; - d,, 
2 1 il 2 1 1 
X5 a d.). Assuming the latter is the smaller of the two, 


a pattern move is made in the direction established by the 


preceding successful moves to the point x! teeanoure’ (3.1). 


1 
At x, the small moves which proved to be successful around 


1) 


x” are repeated and if successful lead to a second pattern 


move to me From xo, small displacements of magnitude d, 


IAAL 











“* 





ae 


FP ajenuge: (Cd ) 


Pattern Search 
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1 coordinate prove to be unsuccessful, whereas 
: : 2 2 2 2 
the displacement d, yields E(x} » X5 + d,) < E(x, 2 x5). 


aalomea the wx 


A pattern move is made therefore along the line which 


1 


passes through the points om = ‘ah 5 


1? + d,) and 


-- 
(x; » X5 d,). 
Note that the step size of the pattern search is made 
progressively larger as the general directions dictated 
by the pattern prove successful. When failure occurs, 


aS at point a mie Step size Of the pattern move is 


a At <a tne  d1spaacenen t 


~e 


reduced as depicted, to point x 


d, from x? along the x, coordinate fail to improve f(x), as 


1 


from x? along the x? COordinate. Reduc- 


2 
fhons in dy and d, are therefore required at this point, 


do displacement d 


and the process starts anew. Ultimately both the dis- 


placement and the pattern step sizes are reduced below 


some preassigned limit, and the search is terminated. 
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